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Abstract. A possible origin of the anomalous dip and bump in the primordial power 
spectrum, which are reconstructed from WMAP data corresponding to the multipole i = 
100 ~ 140 by using the inversion method, is investigated as a consequence of modification of 
scalar field dynamics in the inflation era. Utilizing an analytic formula to handle higher order 
corrections to the slow-roll approximation, we evaluate the relation between a detailed shape 
of infiaton potential and a fine structure in the primordial power spectrum. We conclude 
that it is unlikely to generate the observed dip and bump in the power spectrum by adding 
any features in the infiaton potential. Though we can make a fine enough shape in the power 
spectrum by controlling the feature of the potential, the amplitude of the dip and bump 
becomes too small in that case. 
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1 Introduction 

Inflation has been the focus of attention as a successful scenario which gives not only a natural 
understanding of initial conditions of the standard Big Bang Cosmology, but also a gener- 
ation mechanism of primordial curvature fluctuations, which eventually turn into Cosmic 
Microwave Background (CMB) anisotropics and Large Scale Structure of the universe. Cur- 
rent precise data of cosmological observations, e.g., Wilkinson Microwave Anisotropy Prove 
(WMAP) data [1], predict a nearly scale-invariant power spectrum of primordial curvature 
fluctuations. As is well-known, such kind of the power spectrum can be generated from the 
standard single slow-roll inflation. In keeping with the progress of the cosmological observa- 
tions, however, small characteristic features in the primordial power spectrum, e.g., a lack of 
large-scale power, small dip and bump at multipole £ ~ 20 ~ 40, a fine oscillating structure 
around multipole I ~ 100 ~ 140 and so on[l, 2], have gotten more and more attention. 

In fact, there are many attempts to explain such features in the primordial power 
spectrum by considering a non-trivial dynamics of the scalar field during inflation. For 
example, several authors [4-6] claimed that the small dip and bump at multipole £ w 20 ~ 40 
can be explained due to the temporarily breaking off from the slow-roll dynamics of the 
inflaton field. They employed an inflaton potential with a sharp step proposed by Adams et 



Indeed, the above effective potential can be naturally considered in the two-field inflation 
model with a very small mass difference between two inflaton fields. Here the parameter 
A characterizes the "mass gap". For the small value of A, the general slow-roll approxima- 
tion formula proposed by ref. [11] is applicable to the generation of primordial curvature 
fluctuations. Either employing this approximation or numerically solving the perturbation 
equations, the authors obtained a power spectrum of primordial curvature fluctuations which 
can explain the observed small dip and bump at ^ ~ 20 ~ 40 in the CMB temperature power 



al.[3]: 
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spectrum [4-6]. It should be noted that the width of an oscihation in the power spectrum of 
curvature fluctuations is approximately Alnfc ~ 1 — 1.5 in order to match the WMAP data. 

Recently, a new oscillating structure around multipole i ~ 100 ~ 140 has been reported[7, 
8] . The authors employed the inversion method to reconstruct the primordial curvature fluc- 
tuations from the WMAP data. The width of the structure they found in the primordial 
power spectrum is Aln/c ~ 0.04, which is much finer than the structure at £ ~ 20 ~ 40. 
In order to realize such a fine structure in the primordial power spectrum, Nakashima et 
al.[9] have proposed a sudden change of the sound velocity during inflation induced by a non- 
canonical kinetic term of the inflaton. They found fine oscillations in the primordial power 
spectrum, while these oscillatory features continue on the scales smaller than the Hubble 
radius at the time of the sound velocity change. Therefore such a model may screw up 
the matching between the observational data and theoretical predictions at £ > 200 which 
corresponds to the scale of the first acoustic peak. 

In this paper, we investigate other possibilities of generating the fine structure of the 
primordial power spectrum by modifying the inflaton potential. Employing the general slow- 
roll approximation formula, we can easily find the relation between the fine structure of the 
primordial power spectrum and the evolution of the slow-roll parameters during inflation. 
Here the dynamics of the inflaton determines the evolution of slow-roll parameters. For 
example, it can be shown that the mass transition which is described by eq.(1.2) modifies 
the monotonic evolution of the slow roll parameters and eventually generates a structure in 
the primordial power spectrum. We expect that the shorter the mass transition is, the finer 
the structure becomes. 

This paper is organized as follows. In the next section, we briefly review and investigate 
an analytical formula to obtain the primordial power spectrum, that is, the general slow- 
roll formula. In section 3, based on the formula, we discuss the possibility of generating 
fine structures in the primordial power spectrum due to the modification of the monotonic 
evolution of the slow roll parameters. The last section is devoted to the summary and 
discussion. 



2 General slow-roll formula 



Let us consider a single canonical scalar field (p whose action is given by 

1 



'matter 



(2.1) 



where g^i, is a metric and V{4>) represents the potential of the scalar field. The power 
spectrum of the curvature perturbation TZc induced by quantum fluctuations of the scalar 
field can be defined as 



(7^,(k)7^,(k')) = {27Tf6{k + k')^Vn{k) • 
In the standard single-field slow-roll inflation, the power spectrum is simply given by 



2\ 2 
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where H is the Hubble parameter, (p is the derivative of the inflaton with respect to the cosmic 
time and the subscript aH = k, with a being a scale factor, denotes that the expression is 
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to be evaluated at the time when the scale of interest exits the horizon. The above simple 
formula can be obtained under the assumptions that the slow-roll parameters are small and 
keep almost constant values in time. Even if the case these assumptions are violated, however, 
a general formula of the power spectrum of primordial curvature perturbations proposed by 
Stewart [11] can be applicable, which is refereed as the general slow-role formula. Following 
ref. [11], let us briefly review this formula. 

The time evolution of Tlc{k) can be described by [10] 



where t = J dt/ a is a conformal time, Vk = zTZc{k) and 



+ (k'-z^]vk = 0, (2.4) 
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Using the slow-roll parameters e and rj, we have 
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where 



1 



It should be noted that Eq.(2.6) is applicable for any values of the slow roll parameters, 
while in the standard slow roll approximation, we assume that |e| <^ 1, \rj\ <^ 1 and their 
time derivatives are also small. 

As far as the slow-roll parameters are order of unity or less, we find {l/z){d'^z/dT'^) = 
0{a^H^). Hence, on the super-horizon scales (k <^ aH), the above evolution equation (2.4) 
can be approximately reduced as 

di^Vk 1 d'^z 

From this equation, we can easily find that in the single field inflation the curvature pertur- 
bations TZc, remain constant in time on the super-horizon scales. On the other hand, on the 
sub-horizon scales, namely k » aH, eq. (2.4) can be reduced as 

^ + k'v, ^ . (2.9) 

This equation looks similar to the equation of motion for the harmonic oscillator in the 
Minkowski spacetime and hence the evolution of is decoupled from the background infla- 
tionary dynamics. 

From the above discussion, the power spectrum of TZc strongly depends on the infla- 
tionary dynamics around the horizon crossing time. Therefore, in order to have a specific 
feature in the power spectrum, the modification of slow-roll dynamics has to take place at 
the epoch of horizon crossing for the corresponding scale. 

Now we are at the position to solve the eq.(2.4) by employing the general slow roll 
approximation. Let us first obtain the asymptotic solutions. On the super-horizon scales. 



-3- 



TZc stays constant in time as we have already shown. On the sub-horizon scale, we adopt 
the Bunch-Davies vacuum to provide initial conditions. Accordingly, we have the asymptotic 
solutions as 

1 

-ikT (j^^ _^ —oo) 
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Following ref. [11], we introduce non-dimensional variables y 
also convert z to f{x) as 

f{x) = xz. 



2kvk and x 



Then equation of motion can be rewritten as 



dx'^ 



+ 1 
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where the source function g[x) is defined as 



9{x) = 



dx"^ dx 



(2.13) 



Hence the power spectrum of the comoving curvature perturbation is given by 
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Since the asymptotic behavior x — )• oo of the homogeneous solution of eq. (2.12) has to match 
with the one kr — )• — oo of eq.(2.10), we obtain the homogeneous solution of eq. (2.12) as 
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Therefore, using the Green's function method, we obtain the formal solution of eq. (2.12) 
as 



y{x) = yoix) + 



du 



g{u)y{u) [yl{u)yQ{x) - yl{x)yQ{u)] . 



(2.16) 



Since g{x) can be written in terms of slow-roll parameters and their time derivatives as shown 
below, we can solve the above formal solution iteratively. In fact, under a slow-roll condition; 
e ^ 1, we can rewrite the conformal time as 



k 



1 
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By using this expression, the source function g can be expressed as 



(2.17) 
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Figure 1. The window function We{x). There is a sharp cut-off at x = f and oscillatory behavior is 
shown at x ^ 1. 

Unlike the standard slow-roll approximations, we take into account the time derivatives of 
the parameters, such as e and rj in this general slow-roll formula. 

Then, up to the first order in the source function g, we develop an analytical formula 
of the power spectrum of the curvature perturbations as 



kV- 1 
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(2.19) 



where * denotes the value at the horizon crossing time, = —kr^, = 1, and the window 
function 

We{x) = W{x)-e{x^-x), (2.20) 

with 

= _ 3 cos (2x) _ 3sin(2x) 

2x3 x2 2x ' ^ ' 

^(-) = {i (oSx). (^-^2) 

Figure 1 shows the form of the window function. In this figure, since x = 1 is the horizon 
crossing time for a given wavenumber, x ^ 1 corresponds to the time after the horizon 
exit (on super-horizon scales) and while x ^ 1 to the time before the horizon exit (on sub- 
horizon scales). This function reflects the behavior of the homogeneous solution yo{x). Due 
to the oscillatory behavior of the window function for x ^ 1 and the sharp cut-off for x ^ 1 
coming from the step function, we can easily find that the integral term of the formula mostly 
depends on the behavior of the source function around the horizon cross time. Therefore it is 
clear that the specific feature on the power spectrum induced by the modification of slow-roll 
dynamics appears only at the scales crossing the horizon when the modification takes place. 

If slow-roll conditions are satisfied, the source function g and {df/dx)/f are approxi- 
mately zero. Then, we verify that eq. (2.19) leads to the well-known slow-roll expression (2.3); 
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Figure 2. The source functions g{4>) with all terms 
of the right hand side of eq. (2.13) (red solid line) 
and with only V^^/H"^ (blue dotted line) for the 
bump model. 



Figure 3. The relation between the width of the 
inflaton potential A<pQ and the one of the source 
function A.(f>g for the bump model. 



Therefore, we can conclude that eq. (2.19) is a natural extension of the slow-roll approxi- 
mation to the generalized one. In appendix A, we check the validity of this general slow-roll 
formula. 



3 Can we make fine structures in the primordial power spectrum? 

In the previous section, we have briefly reviewed the general slow-roll formula to evaluate 
the primordial power spectrum in which we iteratively take into account the modification of 
the slow-roll dynamics. The analytical formula eq. (2.19) consists of the window function 
Wg{x) which reflects the behavior of the homogeneous solution yo{x), and the source func- 
tion g{x) related to the inflationary background dynamics through the slow-roll parameters. 
Accordingly, the fine structures of the primordial power spectrum depends on the form of the 
window function and also that of the source function. In this section, using the analytical 
formula, we study the modification of the primordial power spectrum by considering the 
effect of the modification of the slow-roll dynamics. We focus on the relation between the 
fine structures of the primordial power spectrum and the form of the source function g with 
keeping the form of the window function. 

3.1 Source function g 

First, let us briefly mention how the functional form of g can be related with theoretical 
models of inflation. Stewart [11] has calculated the primordial power spectrum by making 
use of the general slow-roll formula with changing the functional form of f{x). In this paper, 
however, we focus on the form of the source function g(x) instead of f{x) which we can 
construct from g by using eq. (2.13). 

In eq. (2.18), we gave the expression for g with respect to the slow-roll parameters defined 
by eq. (2.7). Here, we give another expression for g in terms of derivatives of potential V{(p) 
such as Vdy = dV/dcj) and Vss = cPV/dcj)'^, time derivative of cp and Hubble parameter H as 
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In order to investigate which term of the right hand side of this equation mainly contributes 
to g, in figure 2 we show the source function g and Vfpcp/H'^ for a bump model whose potential 
is given by 

mis = m2 [1 + cexp (-(</. - 0o)VA0g)] . (3.2) 

From this figure, it is found that the form of g is dominated by V^^/H"^. We find that this is 
also true for the case of a step in the potential. Therefore we can conclude that V^^ controls 
g{(f)) when the dynamics of the infiaton modifies due to the fine structure in the potential. 
As is well known, V^^ represents the effective mass of the scalar field. In this sense, it is 
considered that the source function g traces the change of the mass of the infiaton. Of course, 
if the infiaton is coupled with other massive scalar fields during inflation, then the effective 
mass of the infiaton can be affected by the evolution of such massive fields. Hence, such kind 
of fine structures in the primordial power spectrum may provide some clues for high energy 
physics [12-17]. 

Let us now focus on the width of the source function g which may control the width of 
the primordial power spectrum. On the contrary, the width of g should relate with the one of 
the infiaton potential. Figure. 3 shows the relation of the A(/)o to the parameter /\(j)g which is 
defined as the width of g[(f)) as is shown in figure 2. From this figure, it is clear that the width 
of g{(j)) is a monotonic function of the one of the infiaton potential. Hence, when we can 
find the form of the function g through the observations of the primordial power spectrum, 
we can easily reconstruct the potential form of the infiaton. In the following discussion, we 
consider the effect of the anomalous structure of the infiaton potential on the primordial 
power spectrum with changing the functional form of g. 

Before closing this subsection, we would like to remark that we can separate g into 
two components, namely the smooth part (7gsr whose integration in terms of time produces 
a power law component in the primordial power spectrum and the gap part g^^p which 
corresponds to a fine structure in the spectrum as 

g = 5GSR + ffgap- (3.3) 

Note that the function f{x) in eq. (2.19) can be derived by integrating the function g{x) 
from its definition eq. (2.13). Here the power law part can be written by using the spectral 
index Ug, which is defined as V-r. oc as 

^^GSR = ^^-(n.-l). (3.4) 

In this paper, we adopt = 0.96[1]. 

On the other hand, we employ a general form of the infiaton potential with steps and 
bumps to investigate whether we can obtain the observed fine structure at £ ~ 100 ~ 140 
or not in the next subsection. It turns out that the resultant functional form of ^fgap can be 
described by the differential Gaussian function. 

3.2 Fine structures of the power spectrum controlled by g 

In order to investigate the fine structure of the primordial power spectrum due to the mod- 
ification of the slow-roll dynamics, we employ the differential Gaussian function as the gap 
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Figure 4. The source function g for tlie differential Gaussian model with n ~ 1. We set the parameter 
d from 5.0 x 10~^ to 7.5 x 10"^ for the left panel and from 5.0 x 10~^ to 7.5 x 10^"^ for the right panel. 




Figure 5. Ratio between A'P{k) = Vnik) — VcsRik) and VcsRik), where Vn{k) is the full power 
spectrum and T^gsrI^) is the power-law power spectrum which is constructed from gcSR- As d, which 
is the width of g, decreases, the width of the structure in the primordial power spectrum also decreases 
as is shown in the left panel. If c? < 1.0~^, however, the width of the structure no longer decreases. 



part of the source function 5(gap(x), given by 



gga.p{x) = A exp 



d 



[n/2] 

E 



n!(-d)^-" 
i\{n - 2i)\ 



{2{x-h)] 



n—2i 



for n = 1, 2, • • • 



(3.5) 



Here [ra/2] is the Gauss symbol. It should be noticed that n = 1 and n = 2 correspond to 
a step and a bump / dip in the inflaton potential, respectively. Therefore the differential 
Gaussian function for ^(gap is a natural and general extension of the inflaton potential with 
simple fine features. Moreover, the width and height of the features can be easily controlled 
by changing the parameters d and A, respectively. 

Our goal is to explain the newly found fine structure at ^ ~ 100 ~ 140. By employing 
the inversion method, as we have already mentioned, Ichiki et al. [8] found the structure. 
They considered several different features in the primordial power spectrum to match the 
data. They found that a typical width of the fine feature is A In A; = 0.04 and a best-fit 
functional form is v^ one, that is a power spectrum with a peak &t k = and a dip at 
k ^ k^ . 
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Figure 6. The width of the structure in the power 
spectrum in logarithmic space, A In fc as a function 
of the width of the source function d for the n = 1 
case. At d < 1.0 x 10~^, A In A: becomes roughly 
constant. 



Figure 7. The window function (black line) and 
the source function ggap- As the width of the pri- 
mordial power spectrum features does not change, 
the width of the ^gap is less than the that of window 
function. 



In fact, it turns out such a feature in a power spectrum can be realized even if we take 
the simplest functional form, i.e., n = 1 in eq. (3.5) with an appropriate set of parameters. 
However it is not clear whether we can satisfy the height and width of the power spectrum 
simultaneously to fit the observation. Therefore let us begin with the n = 1 case. The source 
function ^gap can be given by 



gga.p{x) = (x - 6)exp 



{x - bf 
d 



(3.6) 



Here, the parameters c and d determine the height and the width of ggap, respectively, and 
b{= — In /co) corresponds to the position of the fine feature on the primordial power spectrum. 
To fit the fine structure at £ ~ 100 ~ 140, we set ko = 0.002 [Mpc~^]. To maximize the 
amplitude of the feature under the assumption 51 ~ 0(1), we set c = 2. 

In figure 4, we show the source function g for several different values of d between 10~^ 
and 10~^. As we decrease the value of d, the width of (/gap becomes finer as we expected. 

In figure 5, we show deviation of the primordial power spectrum V'R.{k) with g = (7gsr + 
ffgap from the power-law power spectrum Pgsr(^) with g = ^gsr for each value of d. Here 
the deviation A'P(A;) = V-jzik) — "PgsrC^)- From this figure, it is found that the form of the 
power spectrum shows a single dip for d > 10^^, while there appears a dip and a bump in 
the power spectrum, that is the v^ type, for d < 10^^. If we push d down to d ~ 0(10"^), 
the power spectrum starts to oscillate. 

The width of the power spectrum A In is defined by the distance between the wave 
numbers where A'P(A;)/Pgsr(^) becomes zero around the maximum value of AV{k)/VGSR{k)- 
In figure 6 , we show the dependence of the width A In A; on the parameter d which controls the 
width of 5gap- In fact, this figure shows that A In k decreases as d decreases for d > 1.0 x 10~^ . 
However, Aln/c starts to saturate for d < 1.0 x 10~^, whose value is about A In A; ~ C(l)- 
It turns out that it is very difficult to realize the observed width AlnA; = 0.04 in this case. 
In order to explain this result, we compare the source function g with the window function 
Wg{x) in figure 7. From this figure, we see that the width of the source function becomes 
smaller than the one of the window function when d < 1.0 x 10~^. Since the width of the 
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Figure 8. Ratio between A7'(fc) = V-jiik) — VcsRik) and VcsRik)- for the n = 2 in the left panel 
and for n = 3 in the right panel. 

primordial power spectrum is determined by the products of g and W0{x), it is mainly deter- 
mined by the width of the window function in this case. This causes the saturation because 
the window function smoothes out tlie finer feature of the source function. On the other 
hand, if d > 1.0 x 10~^ the width of the source function becomes larger than the one of the 
window function. Accordingly, the width of the primordial power spectrum is controlled by 
the one of the source function, which is d. 

Now, let us investigate higher order differential Gaussian models given by eq. (3.5), i.e., 
models with n > 2. In figure 8, deviations of the full primordial power spectrum from the 
power-law one for n = 2 and 3 are shown. We find that the width of the structure becomes 
finer for higher order models, although the saturation for d < 1 x 10~^ can be seen as the 
model with n = 1. In figure 9, the widths of the structure A In A; are plotted as a function of 
d for n = 1, 2 and 3 to see this tendency more clearly. 

We can conclude that, to have a finer structure in the primordial power spectrum, we 
need to employ a higher order differential Gaussian function as a source function. It is 
because a higher order differential Gaussian function has more numbers of oscillation for a 
given d, and these oscillations pick up the larger wave number part of the window function, 
which actually produces finer structure in the power spectrum. 

In figure 10, we push the number of differentiation n much further up to 100. We can 
almost reach the observed value Aln/c = 0.04. However, there is a caveat. If we take n to 
be a large number, the maximum value of AP/'PgsRj which we refere as the amplitude of 
the structure hereafter, becomes smaller. As is shown in the right panel of figure 10, the 
amplitude drops below 10^^ for n = 100, which is way too small to explain the observed 
structure. 

4 Summary and discussion 

In this paper, we have discussed the fine features on the primordial power spectrum which 
can not be characterized by a simple power-law spectrum. Previous works to focus on the 
fine structure of the CMB temperature power spectrum at the multipole / = 20 ~ 40 have 
shown that such an oscillating feature can be realized by the modification of the slow-roll 
dynamics during inflation. Based on this idea, we have examined the possibility of explaining 
a newly found finer feature of the primordial power spectrum which are reconstructed from 
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Figure 9. The width of the structure in the power spectrum in logarithmic space, A In fc as a function 
of the width of the source function d for the n = 2 and 3 cases of differential Gaussian models with 
the n = 1 one (figure 6). It is shown that Alnfc becomes smaller for larger n models. 
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Figure 10. The width of the fine structure Alnfc as the function of the number of differentiation 
?? in the left panel. Even if the n is increased to 100, Alnfc can not reach at the observed value 
Alnfc = 0.04. And the right panel shows that the maximum amplitude of AV{k)/VGSR{k) as the 
function of n. The amplitude drops below 10^^ for n = 100, which is too small to explain the observed 
structure. 



WMAP data corresponding to the multipole i = 100 ~ 140 by using the inversion method 
by ref. [8]. 

Employing the general slow-roll formula presented by ref. [11], we can characterize the 
power spectrum by the window function, which is a oscillating function with a sharp cut off 
corresponding to the horizon crossing scale, and the source function, which is described by 
the slow roll parameters and their time derivatives. For the source function, we have used 
the differential Gaussian function, since it can naturally describe the structures such as step, 
dip and bump in the inflaton potential. Moreover, we can easily control the width of the 
structure and the number of oscillations of the source function. 

We have found that the width of the fine structure Alnfc is controlled by the width 
of the window function in the first place while it also depends on the width of the source 
function d. In fact, since the window function smoothes out the structure of the source 
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function, Alnfc saturates when we decrease the value of d below ~ 1.0 x 10"^. However, we 
have obtained the finer structure by employing the higher order differential Gaussian model 
as the source function. It is because that oscillations of the source function pick up the larger 
wave number part which corresponds to the inner horizon part of the window function. On 
the other hand, the amplitude of the fine structure, which is defined by the maximum value of 
AP(A:)/Pgsr(^); becomes smaller as we increase the index of differentiation n of the source 
function. It turns out that it is hardly possible to recover the observed fine structure at 
i = 100 ~ 140 with simultaneously satisfying the width and the amplitude. 

There might be a possible solution to expalin the fine structure with introducing a higher 
order correction of the general slow-roll formula [18], which allows us to have a large value 
of the source function \g\ > 1. In this case, the source function may dominate the window 
function and control the fine structure of the power spectrum. Such a source function may 
generate a structure with a large amplitude. However, Dvorkin and Hu [5] studied the effect 
of the second order term and found that only a little enhancement of the amplitude takes 
place. 

The only way to explain the observed fine structure might be to have a fine tuned 
source function which perfectly cancels the oscillations of the window function. Of cource, 
such a source function is compleltely ad-hoc. We may conclude that, therefore, it is almost 
impossible to generate the observed fine structure from a simple modification of the slow roll 
dynamics by adding specific features in the inflaton potential. 
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A Validity of general slow-roll formula 

In this appendix, we check the validity of the analytical formula given by eq. (2.19), by 
comparing the approximate solution with a numerical calculation. Here we employ the chaotic 
inflation model with a dip in the potential, whose potential is given by eq. (3.2) with c < 0. 

In the left panels of figure 11, we show the source function g with m = 7.0 x 10~^, 
(t>o = 14.67 and A(f>l = 2.7 x 10"^ for c = -5.6 x 10-^ -5.6 x 10"'' and -1.7 x 10"^ from 
top to bottom. These values of model parameters correspond to the ones employed in ref. 
[4], in which the feature in the WMAP data at ^ = 20 ~ 40 are fitted. Note that we have 
made a conversion of parameters since they assumed a small step in the potential instead 
of a dip with keeping the width and the position of the feature. From these panels, we can 
see that the maximum value of the source function g becomes larger as the absolute value 
of c increases. In the right panels of figure 11, we show the primordial power spectra of 
respective models obtained by making use of the general slow-roll formula (red solid line) 
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and the numerical calculation (blue dotted line) and also the fractional errors between the 
two results. 

From this figure, we find that around the structure of the primordial power spectrum 
the fractional error between the approximate solution and the exact solution obtained by the 
numerical calculation becomes larger as the parameter c increases. If 10% errors are allowed, 
it is confirmed that the general slow-roll approximation is valid for the cases with g = 0{\) 



around the maximum such with c = —5.6 X 10 ^ for the dip model. Therefore, in 

the context we have employed the cases where the maximum value of the source function g 
is order of unity. 
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